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2  MCTEKS 


I.  INTRODUCTION 


This  experiment  demonstrated  the  obtainable  accuracy  using  Doppler  techniques  for  estimat¬ 
ing  position  changes  of  a  Global  Positioning  System  receiver  antenna. 

Accurate  Doppler  measurements  were  taken  using  the  NAVSTAR  Geodetic  Receiver  System 
(NGRS),  developed  at  the  Naval  Surface  Weapons  Center.  The  system  uses  a  GPS  receiver  built  by 
Stanford  Telecommunications  Incorporated  of  Sunnyvale,  California,  and  a  Hewlett-Packard 
506/A  high  performance  cesium  clock. 

The  experiment  consisted  of  periodically  moving  the  receiver  antenna  from  its  original  position 
(location  0  of  Figure  1 )  to  approximately  2  meters  (m)  east  of  the  original  position  (location 
1)  and  then  to  a  third  position  approximately  2  m  south  of  the  original  position  (location  2). 
Figure  1  presents  the  surveyed  relative  locations. 
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Figure  1.  Relative  Antenna  Locations 
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Each  antenna  move  was  timed  to  be  within  a  one-minute  Doppler  count  interval.  The  moves 
occurred  every  fifteenth  interval,  except  during  a  period  between  tracking  of  the  two  satellites 
used  to  collect  data.  The  moves  were  repeated  until  at  least  six  measurements  were  made  to  each 
location.  These  measurements  were  used  to  estimate  the  three-dimensional  position  change  of  the 
antenna.  This  emulated  a  multichannel  receiver  tracking  six  satellites  simultaneously. 

The  next  section  presents  the  calculated  change  in  range  due  to  the  antenna  position  change 
only.  Section  III  gives  the  procedure  for  obtaining  the  observed  change  in  range  and  compares  it 
with  the  calculated  values.  Section  IV  describes  the  procedure  for  estimating  the  three-dimensional 
antenna  position  change,  and  Section  V  presents  the  results.  The  accuracies  of  the  experiment 
are  discussed  in  the  final  section. 


U.  CALCULATED  CHANGE  IN  RANGE 


The  calculated  change  in  range  from  the  satellite  to  the  antenna  due  to  the  antenna  position 
change  AR^  is  shown  geometrically  in  Figure  2.  Parallel  lines  are  assumed  from  both  locations  to 
the  satellite  since  the  distance  to  the  satellite  is  much  larger  than  the  distance  moved.  The  figure 
shows  the  Ad  baseline  direction  to  be  a  degrees  from  north.  The  values  of  a  for  the  three  position 
changes  of  the  experiment  are  given  in  Table  1.  For  example,  for  location  0  to  location  1  ct  is  90°. 
The  survey  values  for  the  distance  moved  in  the  a  direction.  Ad,  and  the  perpendicular  distance, 
Ap,  are  also  given  in  the  table.  The  antenna  was  moved  in  the  horizontal  plane  so  that  the  change 
in  altitude.  Ah,  is  zero.  The  accuracy  of  the  survey  values  should  be  within  1 .0  centimeter  (cm). 


Figure  2.  Geometric  Representation  of  Position  Change 
from  Location  0  to  Location  1 


Table  1.  Surveyed  Change  in  Antenna  Position  (m) 


Location 

Change 

Of 

(°) 

Ad* 

Ap* 

Ah* 

Oto  1 

90 

2.000 

0.133 

.000 

1  to  2 

225 

2.779 

0.238 

.000 

2  to  0 

0 

2.000 

0.203 

.000 

Accuracy  to  within  1  cm 


The  relationship  for  the  change  in  range  is  found  from  a  first-order  expansion  about  the  origin 
(location  0). 


AR^(t.)=  (t.,d  + Ad,  p  + Ap,  h-fAh)-  (t,d,  p,  h)  (1) 

3Rc  3Rc  aRc 

(tj.d,  p,  h)Ad+- — (t.,  d,  p,  h)  Ap +  — (t.,  d,  p,  h)  Ah.  (2) 

od  op  '  oh 


From  Figure  2,  the  change  in  calculated  range  for  small  position  changes  from  the  origin  are 
given  as  follows 


aR^itj,  d,  p,  h) 

8d  Id  =  p  =  h  =  0 


cos  0(tj)  cos  E(tj) 


(3) 


BRf  (tj,  d.  p,  h) 


ap 


W=p=h=0 


sin  0(t)  cos  E(t.) 


(4) 


aR(.  (tj.  d.  p,  h) 

ah  d=p=h=0 


=  -  sin  E(t) 


(5) 


i 

! 
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where 


t  =  time  of  the  expansion 
l'(tj)  =  elevation  to  the  satellite 
A^ttj)  =  azimuth  to  the  satellite 
a  =  direction  of  the  base  line 
0(tj)  =  a  -  A^(f) 

Ad  =  survey  value  for  change  in  position  in  the  q  direction 
Ap  =  survey  value  for  change  in  position  in  the  ck  -  90°  direction 
Ah  =  survey  value  for  change  in  position  in  the  vertical  direction 


III.  OBSERVED  CHANGE  IN  RANGE 


I'here  are  several  steps  in  the  procedure  to  compute  the  oLserved  change  in  range.  AR^,  The 
Doppler  measurement  is  obtained  by  integrating  over  a  lime  interval,  ftj.  t^J,  the  received  fre¬ 
quency,  f^,  subtracted  from  a  precise  ground  station  frequency,  f^.  The  received  frequency  is 
the  sum  of  the  transmitted  frequency,  f.j,,  and  the  Doppler  effect.  Therefore,  the  Doppler  measure¬ 
ment  is 


AN  = 


t 

f 


t 

I 


tj  R(t) 


dt 


(bl 


=  (f^  -  -  t,  )  +  ^  [R(t2)-  R(t,  ))  (7) 

where  c  is  the  speed  of  light  and  R(t)  is  the  instantaneous  range  from  the  satellite  to  the  station 
antenna.  Note,  the  Doppler  measurement  is  a  difference  in  ranges  at  the  beginning  and  end  times 
('t  the  Dopplv’r  count  interval.  The  range  at  any  time  between  these  times  is  immaterial.  For  this 
reason,  the  station  antenna  for  the  experiment  was  in  place  at  its  first  location  at  t| ,  moved  after 
tj  .  and  pul  in  place  at  its  second  location  before  t^ .  The  start  and  end  times  for  each  position 
cliange,  as  well  as  their  corresponding  azimuth  and  elevations,  are  given  in  Table  2. 

file  (1P.S  satellites  transmit  two  signals  at  frequencies  that  enable  a  correction  to  be  made  for 
the  first-order  ionospheric  effect.  After  ionospheric  correction  of  the  range  difference  biased 
observed  ranges  Rp(tj)  are  formed  by  adding  successive  range  differences  to  a  nominal  range 
bias.  Next,  trajectories  are  computed  for  the  satellites,  and  calculated  ranges  R^(tj,  d.  p,  h)  are 
determined  to  location  0  of  Figure  1.  Next,  a  least-squares  fit  is  made  using  the  observed  and  the 
calculated  ranges  to  determine  the  unknown  bias  Rg  and  improvements  to  six  orbit  parameters 
and  tropospheric  refraction  corrections.  The  residuals  of  the  fit  are  plotted  in  Figure  3  for  satellite 
.“i  and  in  Figure  4  for  satellite  6. 
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Table  2.  Azimuths  and  Elevations 


Satellite 

Location  Number 

‘i-i* 

Azimuth 

Elevation 

E(ti.,)(°) 

‘i* 

Azimuth 

A,(tj)(°) 

0-1  5 

40500 

342.8 

54.4 

40560 

343.6 

43200 

23.2 

60.6 

43260 

24.2 

45900 

70.2 

59.1 

45960 

71.2 

48600 

107.2 

48.2 

48660 

107.9 

51300 

130.5 

31.9 

51360 

130.9 

6 

60300 

61.3 

23.0 

60360 

60.8 

63000 

39.7 

19.7 

63060 

39.3 

1-2  5 

41400 

355.1 

57.0 

41460 

355.9 

44100 

38.8 

61.2 

44160 

39.9 

46800 

84.2 

56.4 

46860 

85.1 

49500 

1 16.3 

43.2 

49560 

1 16.8 

52200 

136.1 

25.9 

52260 

136.5 

6 

61200 

76.2 

19.5 

61260 

75.7 

63900 

53.7 

23.1 

63960 

53.2 

2-0  5 

42300 

8.5 

59.1 

42360 

9.4 

45000 

54.8 

60.7 

45060 

55.8 

47700 

96.6 

52.7 

47760 

97.3 

50400 

124.0 

37.7 

50460 

124.5 

6 

62100 

68.9 

21.8 

62160 

68.4 

64800 

46.4 

21.9 

64860 

45.9 

•Year  1980,  Day  1 14  (GMT);  unit,  second  of  week 


Elevation 


54.6 

60.7 

58.9 

47.9 

31.5 

23.1 

19.5 

57. 

61. 

56.1 

42.8 

25.5 

19.6 
23.0 

59.3 

60.6 

52.4 
37.3 

21.9 
21.8 
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Figure  3.  Range  ResiduaU  Versus  Time  for  Satellite  5 
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Range  Residuals  Versus  Time  for  Satellite  5  (Continued) 
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These  plots  consist  of  the  least-square  range  residuals  in  meters  on  the  abscissa  and  sequential 
time  on  the  ordinate,  increasing  downward.  Since  the  computational  fit  procedure  assumed  that 
the  antenna  did  not  move,  discontinuties  occur  at  the  time  intervals  when  the  antenna  is  moved 
from  each  location.  The  location  changes  are  given  in  the  right-hand  column  of  the  figures.  These 
discontinuties  are  the  changes  in  range  due  to  the  antenna  position  change  plus  measurement 
errors.  While  tracking  satellite  5,  the  antenna  was  initially  placed  at  location  0  of  Figure  1 .  The 
first  discontinuty  of  Figure  3  represents  the  first  move  of  the  antenna  from  location  0  to  location  1 . 
The  antenna  remained  at  location  1  for  fifteen  1 -minute  Doppler  count  intervals  and  was  then 
moved  to  location  2.  These  moves  were  continued  clockwise  around  Figure  1  and  are  denoted  in 
Figure  3  until  the  end  of  the  pass.  The  position  changes  were  then  continued  for  satellite  6  for  its 
entire  pass,  as  shown  in  Figure  4. 


The  observed  change  in  range  is  taken  as  the  difference  in  the  residuals  before  and  after  the 
position  change;  i.e. 

'I'  '’■"'L.,.,.,  -  "‘o'Vi*-  ^  Oi-I.l-  P. 

(8) 


where  the  antenna  move  occurred  after  tj_|  and  before  tj.  The  observed  and  calculated  change 
in  range  values  are  given  both  in  Figures  3  and  4  at  the  discontinuities  and  in  Table  3.  In  Figures  3 
and  4  the  observed  values  are  above  the  discontinuities  and  the  calculated  values  are  below,  in 
parenthesis.  Both  values  are  in  units  of  meters.  Also  given  in  Table  3  are  the  satellite  numbers  and 
the  start  time  of  the  Doppler  count  interval.  The  observed  and  calculated  changes  in  range  are 
differenced  for  each  measurement.  The  mean  and  standard  deviations  of  (AR^  -  AR^)  for  each 
location  change  are  also  given  in  Table  3.  The  combined  mean  is  0.0016  meters  and  the  combined 
standard  deviation  is  0.0573  meters. 


Note  that  the  standard  deviation  of  the  location  change  2-0  was  significantly  smaller  than  the 
two  other  location  changes.  By  observing  Figures  3  and  4  note  that  the  location  changes  2-0 
always  occured  at  least  one-half  hour  from  either  end  of  the  satellite  pass.  The  processing 
algorithm,  from  which  these  residuals  were  obtained,  attempts  to  produce  a  zero  mean  result. 
Consequently,  there  is  a  steep  slope  to  the  residuals  at  the  begining  and  end  of  the  passes.  This 
phenomenon  is  very  severe  for  the  short  pass  of  satellite  6,  less  so  for  satellite  5.  but  casts  some 
doubt  on  the  accuracy  of  the  first  and  last  observations  in  both  passes. 

Ihe  individual  observed  minus  calculated  range  differences  are  divided  by  the  cosine  of  the 
corresponding  elevation  angle  and  plotted  along  the  azimith  to  the  satellite  from  location  I  (see 
Figure  5).  These  represent  a  zero  elevation  plot  of  the  range  difference  errors.  If  the  survey  and 
measurements  were  perfect,  all  the  difference  values  would  be  at  the  antenna  location.  This  plot 
depicts  the  orientation  of  the  satellites  with  respect  to  the  base  line  between  location  0  and  1 . 
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Table  3.  Observed  and  Calculated  Change  in  Range 


Satellite 


location  Number 

‘i-t* 

ARflltj,  t,_, 

0-1  6 

60300 

-1.692 

63000 

-1.207 

5 

40500 

0.223 

43200 

-  .454 

45900 

-  .997 

48600 

-1.168 

51300 

-1  1 10 

1-2  6 

61200 

1.980 

63900 

2.490 

5 

41400 

1.222 

44100 

1.299 

46800 

0.997 

49500 

0.474 

52200 

-  .227 

2-0  6 

62100 

-  .480 

64800 

-1.1 10 

5 

42300 

-1.004 

45000 

-  .451 

47700 

.240 

50400 

1.030 

ARj(m) 

AR„-ARj(m) 

-1 .674 

-.018 

-1.299 

0.092 

0.270 

-.047 

ARo- 

AR^ 

=  0.028 

-  .447 

-.007 

-  .989 

-.008 

®ARn 

-AR 

=  0.062 

-1.247 

0.079 

0 

c 

-1.218 

0.108 

2.125 

-.145 

2.494 

-.004 

1.074 

0.048 

ARq- 

AR 

=  -.033 

1.333 

-.044 

C 

1.109 

-.112 

R 

-AR 

=  .072 

0.485 

0.011 

c 

c 

-  .262 

0.035 

-  .493 

0.013 

-1.143 

0.033 

AR-- 

■AR^ 

=  .01 1 

-1 .000 

-.004 

V 

-  .483 

0.032 

-A  R 

=  .021 

0.262 

-.022 

0 

c 

1.018 

0.012 

•Yejf  1980,  Day  1 14  (GMT);  unit  second  of  week. 
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Figure  5.  Change  in  Range  Residuals  Along  Azimuth  to  the  Satellite  with  Elevation  Effect  Removed 
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IV.  ESTIMATION  OF  CHANGE  IN  POSITION 

This  section  describes  the  procedures  for  estimating  the  change  in  antenna  position  using  the 
observed  change  in  ranges,  described  in  the  previous  section,  as  data.  The  relationship  for 
ARgCtj,  tj  , )  of  (8)  is  used  to  create  a  linear  equation  in  terms  of  the  unknown  position  changes 
Ad,  Ap,  and  Ah.  This  is  done  by  substituting  the  appropriate  relationship  of  either 
d  +  Ad,  p  +  Ap,  h  +  Ah)  or  R(,(tj,  d  +  Ad,  p  +  Ap,  h  +  Ah),  evaluated  at  d  =  p  =  h  =  0,  into  (8) 
to  match  the  specific  location  change.  This  is  explained  below  for  the  three  location  changes. 

A.  LOCATIONS  0  TO  1 

From  (1 )  and  (2),  with  evaluations  at  the  origin  or  location  0, 

R(.(tj,  d,  p,  h)  =  Rg(tp  d  +  Ad,  p  +  Ap,  h  +  Ah) 

3R-  aRe  aR- 

-  —  (t.,d.  p,  h)-— (t  d,  p,h)-—  (t,  d,p,  h)  (9) 

.  ad  '  ap  '  ah  * 

is  substituted  into  (8).  Using  the  partial  derivatives  of  (3),  (4),  and  (5)  this  reduces  to  the 
following  linear  equation  in  terms  of  position  change  parameters  Ad^^j ,  Ap^j ,  and  Ah^  j ,  which 
are  to  be  estimated.  Here,  the  subscripts  denote  the  specific  location  change. 

ARo(t.,  )  =  -  cos  6  (t.)  cos  E(t. )  Ad^,,  -  sin  0  (t^)  cos  E(t.)ApQ, 

-  sin  K  (tj)  Ah„,  +e(tj.  tj_,  >  (10) 

where  the  measurement  error  is 

e( t. ,  t,.  ,  )  =  R<;  '( t, )  -  R|,' '  (tj,  d  +  Ad„ , ,  p  +  App , .  h  +  Ahp , ) 

-  |Rj"’(f_,  )-  Rj,'»»(tj.,.d,p.h)l  (11) 

and  where  the  superscripts  indicate  the  location  of  the  evaluation. 

B.  LOCATIONS  2  TO  0 

For  this  case  the  antenna  is  at  location  2  at  t._j .  Therefore,  the  expansion  of  (1)  and 
(2)  is  evaluated  at  t^  , .  As  above  for  (9), 

R^(t.  I ,  d,  p,  h )  =  R^(tj_  I .  d  +  ^d,  p  +  Ap,  h  +  Ah) 


II 


(12) 


3R, 


dR 


3R, 


3p 


dh 


The  following  equations  are  similar  to  equations  ( 1 0)  and  (11)  except  for  parameters  q  ,  Ap^  q 
and  Ahjfj. 

^Rq  (‘i-  tj_,)  =  cosd  (t._j )  cos  E(tj_, )  Adjo  +  sin  0  (tj_j)cos  E(t._j)  Ap^^, 


+  sin  E(tj_,)Ah2(, +e(t.,tj_,)  (13) 

where  the  measurement  error  is 

e(VVi)  =  R^0)(tj)  -  R<«)(tj.d,  p.h) 

-  ( R(2>  (tj_, )  -  R^)  (V, .  d  +  Ad,„,  p  +  Ap^o.  h  +  Ahjo)]  (14) 

and  where  the  superscripts  indicate  the  location  of  the  evaluation. 

C.  LOCATIONS  I  TO  2 

For  this  case  both  Equations  (9)  and  (12)  are  substituted  into  (8)  since  both  locations 
are  away  from  the  origin.  This  produces  linear  equations  in  terms  of  the  position  parameters  to 
be  estimated  Adj  AP|  j  and  Ahj  ^  and  in  terms  of  Ad^, ,  Ap^j  and  Ah^, .  Survey  values  are 
used  for  Ad^j ,  Ap^j  and  Ah^j  since  the  antenna  positions  are  assumed  to  be  known  perfectly 
before  each  move.  The  resulting  relationship  is 

ARo(tj.  tj.,  )  =  -  cos  e  (t.)cos  E(t.)  Ad,2  -  sin  9  (t.)  E(tj)  Apj^ 

-  sin  E(tj)  Ahj  2  +  (cos  0  (t_, )  cos  E(t._, ) 

-  SOS  9  (tj)  cos  E(t.)l  Adgi  +  (sin  0(tj_, )  cos  E(tj_, ) 

-sin  0(tj)  cos  E(tj)l  Apg  j  +  [sin  E(t._j ) 

-  sin  E(tj))  Ahg,  +e(t,tj_,)  (15) 

where  the  measurement  error  is 

e(ti,t._,)=R^2)(t,)_  R(2)(t.^d  +  Ado,  +Ad,2,P  +  Apo,  +Ap,2, 

h  +  Ah„,  +Ah,2)-(R<'>(ti.,)-R<'>(V,,d  +  Ad„,, 
p  +  Apo,,h  +  Ah,„)J 

and  again  where  the  superscripts  indicate  the  location  of  the  evaluation. 
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(16) 


V.  RESULTS 


The  linear  equations  of  (10),  (13),  and  (15)  were  used  to  estimate  the  three-dimensional 
position  changes  in  terms  of  Ad,  Ap,  and  Ah,  which  are  defined  in  Section  2. 

Observed  changes  in  range  RQ(tj,  tj_j)  of  Table  3  and  appropriate  azimuths  and  elevations 
of  Table  2  were  substituted  into  these  equations  to  form  an  over-determined  set  of  linear 
equations.  Six  measurements  were  used  for  each  position  change  in  order  to  emulate  a  multi- 
channeled  receiver  tracking  six  satellites.  (For  locations  2  to  0  all  six  measurements  were  used; 
for  both  locations  0  to  1  and  1  to  2,  which  had  seven  measurements,  one  point  near  the  end  of  the 
pass  was  eliminated.)  The  least-squares  solutions  for  the  three  location  changes  are  given  in  Table  4 
together  with  the  corresponding  surveyed  values  of  the  position  changes.  The  estimation  errors, 
the  differences  between  the  estimated  and  suveyed  values,  are  also  given.  The  magnitudes  of  these 
errors  are  3.0,  5.4,  amd  6.3  cm,  for  position  changes  2  to  0,  1  to  2,  and  0  to  1,  respectively.  This 
demonstrates  the  subdecimeter  accuracy  of  the  receiver. 

The  standard  deviations  on  the  individual  estimates  are  also  presented  in  Table  4.  These  were 
determined  (see  Reference  1)  by  weighting  the  square  roots  of  the  diagonal  terms  of  the  3  x  3 
inverse  of  the  parameter  coefficient  matrix  by  standard  deviations  of  (ARq  -  AR^),  given  in 
Table  4  for  each  location  change.  The  estimates  are  consistent  with  their  standard  deviations 
except  that  the  altitude  errors  are  smaller,  compared  to  the  Ad  and  Ap  errors,  than  the  standard 
deviations  indicate. 

The  last  three  columns  of  Table  4  have  the  estimation  errors  for  Ad  and  Ap  converted  to 
errors  in  north  and  east  directions  and  the  altitude  error  repeated.  The  estimation  error  for  “com¬ 
pleting  the  loop”  is  the  sum  of  direction  errors  for  each  leg.  Here,  the  antenna  location  is  assumed 
to  be  known  perfectly  at  the  beginning  of  each  leg.  The  resulting  error  is  3  cm  north,  -  4  cm 
south  and  less  than  -  1  cm  in  altitude. 
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VI.  CONCLUSION 


The  experiment  demonstrated  subdecimeter  change  in  range  measurement  capability  and 
three-dimensional  antenna  position  change  estimation  accuracy.  Note  that  very  few  errors,  which 
would  exist  for  an  absolute  position  test  or  a  two-receiver  relative  position  test  (see  Reference  2), 
entered  into  this  change  in  position  experiment.  Due  to  the  small  baseline  absolute  orbit  errors, 
atmospheric  errors  and  clock  bias  errors  have  very  little  effect  on  the  results.  The  three  main 
contributing  errors  of  the  experiment  are  jitter  in  the  measurement  of  number  of  cycles  in  the 
Doppler  count,  drift  in  the  satellite  and  receiver  clocks  over  the  1-minute  Doppler  count  interval, 
and  change  in  orbit  error  over  the  Doppler  count  interval.  The  standard  deviation  of  the  Jitter  in 
the  measurement  and  the  combined  clock  drifts  were  expected  to  be  about  1  and  2.5  cm, 
respectively.  The  change  in  orbit  error  appeared  to  be  small  over  the  major  portion  of  the 
trajectory,  but  due  to  the  processing  anomoly  discussed  in  Section  III,  the  residuals  may  not 
be  as  accurate  at  the  pass  extremes  as  in  the  center.  Note  that  for  the  2  to  0  location  change, 
measurement  data  was  not  taken  at  the  beginning  or  end  of  the  pass.  Consequently,  the  accuracy 
of  the  2  to  0  location  change  results  are  significantly  better  for  this  change  in  location. 

The  sampled  standard  deviations  of  the  change  in  range  measurement  errors  for  locations  2  to 
0,  1  to  2  and  1  to  0  were  2.1,  7.2,  and  6.3  cm,  respectively.  These  standard  deviations  are  based 
only  on  six  or  seven  measurements.  The  combined  mean  and  standard  deviation  based  on  20 
measurements  is  0.2  and  5.7  cm,  respectively. 

The  three-dimensional  position  change  estimation  error  is  the  difference  between  the  estimated 
and  surveyed  values.  The  magnitudes  of  the  estimation  error  for  position  changes  2  to  0.  1  to  2. 
and  0  to  1  were  3.0,  5.4,  and  6.3  cna,  respectively. 

These  are  fairly  impressive  results  even  though  the  baseline  was  very  small  and  few  errors 
entered  into  the  experiment.  Further  testing  will  be  done  for  larger  baselines  and  with  more  than 
one  receiver. 
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